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Abstract 

We present the QCD corrections to Rb and to the AB — 1 effective Hamihonian in models with 
a second Higgs field that couples to the quarks respecting the criterion of Minimal Flavor Violation, 
thus belonging either to the (l,2)i/2 or to the (8, 2)i/2 representation of SU{3) x SU{2) x U{1). 
After the inclusion of the QCD corrections, the prediction for Rt becomes practically insensitive 
to the choice of renormalization scheme for the top mass, which for the type-I and type-II models 
translates in a more robust lower bound on tan /?. The QCD-corrected determinations of Rb and 
BR{B — > Xs 7) are used to discuss the constraints on the couplings of a (colored) charged Higgs 
boson to top and bottom quarks. 



e-mail : 

degrassi(§f is .uniromaS . it 
slavichOlpthe . jussieu. fr 



1 Introduction 



One of the main goals of the present experimental program at the Tevatron and at the Large Hadron 
Collider (LHC) is the search for the Higgs boson(s) in order to elucidate the mechanism of electroweak 
symmetry breaking (EWSB). In the Standard Model (SM), the latter is realized in the most economical 
way via a single Higgs doublet. This minimal realization predicts a single neutral Higgs boson, whose 
mass can be constrained, from electroweak precision data and the direct search limit from LEP, to 
be lighter than ~ 200 GeV. However, at the moment, there is no direct experimental evidence for a 
neutral Higgs boson or any other scalar particle like, for example, a charged boson that can be present 
in models with a nonminimal Higgs sector. 

The LHC is going to explore physics up to the TeV scale in order to search for the Higgs boson, 
as well as for any new phenomenon that would confirm the widespread expectation that the picture 
of particle physics in terms of the SM is incomplete. However, new particles with mass in the TeV 
range that couple to quarks at the tree level can modify the predictions for Flavor Changing Neutral 
Current (FCNC) processes. Thus any extension of the SM, starting from the simplest we can think of, 
a Two-Higgs-Doublet Model (2HDM), needs to face the problem of avoiding conflicts with the strict 
limits on FCNC processes. 

Glashow and Weinberg addressed this issue proposing the principle of Natural Flavor Conservation 
(NFC) [1], which requires that the matrices of Yukawa couplings to up and down quarks for all the 
Higgs fields be diagonal in the basis where the quark mass matrices, M^'^, are diagonal. This implies 
that, with the exception of models with vectorlike quarks that mix with the ordinary ones, NFC 
models do not have tree-level FCNC couplings. In the 2HDM case, NFC can be realized imposing the 
sufficient condition that each of the quark mass matrices is obtained from a single Higgs field. This 
can be enforced via a symmetry that acts differently on the two Higgs doublets, leading to two 
possibilities usually referred to as type-I (i.e., the model in which both up and down quarks get their 
masses from Yukawa couplings to the same Higgs doublet) and type-II models (where up and down 
quarks get their masses from Yukawa couplings to different Higgs doublets). 

A less-restrictive way to suppress FCNC processes, still avoiding conflict with the experimental 
bounds, is to consider the criterion of Minimal Flavor Violation (MFV) [2|, which amounts to assuming 
that all the new flavor-changing transitions, including those mediated at the tree level by electrically 
neutral particles, are controlled by the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Thus, the MFV 
hypothesis requires that all the flavor-violating interactions of the new particles present at the TeV 
scale be linked to the known structure of the Yukawa couplings. 

The enforcement of the MFV hypothesis to the case of multi-Higgs models has been recently 
investigated by several groups [3llll|5j. In particular, in ref. [3] it has been shown, via group-theoretic 
arguments, that the MFV hypothesis can be enforced requiring that all the Higgs Yukawa-coupling 
matrices be composed from the pair of matrices and that are responsible for the breaking of 
the SU{3)q^ X SU{3)uji x SU{3)uj^ quark flavor symmetry. This requirement restricts the allowed 
SU{3) X SU (2) X U (1) representations of the Higgs flelds that can couple to the quarks to either be 
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equal to that of the SM Higgs field, i.e. (l,2)x/2) or transform as (8,2)i/2- Examples of the former 
case, besides the NFC type I and II models, are the aligned model of ref. [5] or the class of 2IIDM 
presented in ref. [6]. The latter case is quite different, because the second field does not acquire a 
vacuum expectation value (vev) and does not mix with the SM Higgs field. Thus the scalar spectrum 
of this model contains a CP-even, color-singlet Higgs boson (the usual SM one) and three color-octet 
particles, one CP-even, one CP-odd and one electrically charged, which are split in mass proportionally 
to the SM-Higgs vev [3]. These colored scalar particles give rise to an interesting phenomenology for 
the LHC, not only because - if they are not too heavy - they can be directly produced, but also 
because their indirect effects can influence flavor, electroweak and Higgs physics [3l Fft [8l [9| [TO]. 

Models with a second Higgs doublet present a new and interesting phenomenology, in particular 
related to the presence of a charged scalar. In the flavor sector, decays mediated by a weak charged 
current are the natural place where effects due to a charged Higgs boson, H^, can show up. In the 
electroweak sector the observable Rf, = T{Z — )• bb)/T{Z — )• hadrons) shows a sensitivity to because 
of the specific vertex corrections introduced by the interaction of with the top and bottom quarks. 
Many studies (for the most recent see, e.g., refs. [TTl [121 [13] ) used various combinations of fiavor and 
electroweak observables to constrain the parameter space of the type-II 2HDM, which garnered most 
of the attention because of its property of having the same Higgs-sector realization as the Minimal 
Supersymmetric Standard Model (MSSM). Other studies [HI |T5l [16] explored the parameter space 
of 2HDMs unconstrained by a Z2 symmetry, with the second Higgs doublet still in the (1, 2)1/2 
representation. 

The theoretical accuracy of the predictions in the 2HDM with MFV is not yet at the same level 
as in the SM. Here we take a first step in improving this situation, by (re) considering the QCD 
corrections to two observables, Rb and BR(i? — )• ^^7), which allow to set important constraints 
on the mass and couplings of the charged scalar. The QCD corrections can play a relevant role in 
reducing the error of the theoretical predictions, a well-known example of this fact being indeed the 
radiative decay of the B meson. The present knowledge in the 2HDMs of the two observables we are 
considering can be summarized in this way: in the case of models with a second Higgs doublet in the 
(1,2)1/2 representation, the complete one-loop calculation of T{Z — )• bb) is available [T7], but (to our 
knowledge) no QCD corrections to the charged-scalar contributions are known. The process i? — )• Xj 7 
is instead fully known at the Next-to-Leading Order (NLO) in QCD [181 [H [201 [21] . The case with 
colored scalars in the adjoint representation of SU{3) is less studied. The one-loop charged-scalar 
contribution to T{Z — t- bb) was reported in ref. [7| (see also ref. [22]) while for the radiative decay of 
the B meson only a partial result for the Leading Order (LO) Wilson coefficients of the magnetic and 
chromo-magnetic operators has been presented [3]. 

In this paper we present the QCD corrections to the contribution to r(Z — )• 66) of a charged scalar 
in either the (1, 2)^/2 or the (8,2)1/2 representation. Concerning B ^ Xg^, we compute the 0{as) 
contribution to the Wilson coefficients due to a colored charged scalar in the (8, 2)1/2 representation. 
This is the missing piece to achieve NLO predictions for BR(i? — > Xs 7) for all 2HDMs with MFV. 
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Because of the specific interactions of the colored scalar with the gluons, the Wilson coefficients cannot 
be simply obtained by an appropriate color-factor rescaling of the known (1,2)^/2 result. 

The paper is organized as follows: in the next section we discuss the couplings of the charged Higgs 
boson in the different realizations of the 2HDM with MFV. In section 3 we present the results for the 
QCD-corrected contribution to T{Z — )• bb) due to a charged scalar, covering both cases of color-singlet 
and color-octet particle. We show that, after the inclusion of the QCD correction, the prediction for 
Rh is practically insensitive to choice of an MS or on-shell (OS) renormalization scheme for the top 
mass. The bounds set by Rb on the tbH^ coupling are also shown. Section 4 contains the result for 
the NLO Wilson coefficients in the AB = 1 effective Hamiltonian (the explicit analytic expressions are 
presented in the appendix). The known results for the colorless 2HDMs are recovered, while the case 
of a colored charged scalar is fully new. The restrictions imposed by S —t- X, 7 on the charged-scalar 
interaction with quarks are discussed. Finally, in section 5 we present our conclusions. 

2 Minimally flavor violating 2HDMs 

In a generic 2HDM it is always possible to rotate the two Higgs fields to a basis in which only one of 
them, which we denote as gets a vev [23]. In this basis, we write the Yukawa interactions of the 
Higgs fields with the quarks as 

-£y = QL^iY^UR + QL^iY^'dR + QL^^'^T^^^Y^UR + qL'^i^'^T^R^Y^'dii + h.c. , (1) 

where = ia2^*, and the Yukawa couplings y^'^ are 3x3 matrices in fiavor space such that 
M^'^ = y^'^ i^i)- The possibility of a colored second Higgs doublet is encoded in the matrices 
T^^^ that act on the quark fields. In the usual colorless 2HDM Tr is equal to the identity matrix in 
color space. On the other hand, for a colored Higgs doublet in the adjoint representation of SU (3) 
= Tp (a = 1,8), the matrices of the fundamental representation. The MFV condition amounts to 
requiring that the Yukawa coupling matrices of the second doublet, Y^'^, be composed of combinations 
of the matrices y^'^, and transform under the SU{3)q^ x SU{3)uj^ x SU{2>)ijjj quark flavor symmetry 
in the same way as Y^'^ themselves. We can decompose the matrices Y^'^ as 

Y^ = (1 + 6„y^y^t + . . .) yf/ , Y^ = Aa{i + 6.y^y^t + . . .) y^ , (2) 

where in principle A^^d and are arbitrary complex coefficients. The ellipses in eq. ([2]) denote terms 
involving powers of y^y^t ^s well as terms involving higher powers of Y^Y^"^ . In the following, 
we will assume that the only significant deviations from proportionality between Y^'^ and y^'^ are 
controlled by the Yukawa coupling of the top quark, and that terms involving higher powers of the 
Yukawa matrices are suppressed (e.g., because they are generated at higher loops). If we further 
require that there are no new sources of CP violation apart from the complex phase in the CKM 
matrix, the coefficients ^n,d and Cu^d must be real. Finally, the case = = corresponds to the 
NFC situation in which the Yukawa matrices of both Higgs doublets are aligned in flavor space. 
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The processes Z ^ bb and B ^ Xgj that we will consider in sections 3 and 4 involve loops with 
a charged Higgs boson and a top quark. Under the assumptions implicit in eq. ([2]), the interaction 
between the quarks and is controlled by the Lagrangian 

= -71^ E '^R "^-^ - md, ^-^) + h.c. , (3) 

where g is the SU{2) coupling constant, i,j are generation indices, m^.d quark masses, V is the 
CKM matrix. The family-dependent couplings ^ read 

Ai,d = ^u,d (^1 + eu,d ^ fe) , (4) 

where v = ($1). It appears from eq. (j4]) that, when we neglect the masses of the light quarks, the 
effect on the charged-Higgs couplings arising from the y^y^t terms in eq. ([2]) is limited to a shift in 
the couplings involving the top quark. Since those are the only couplings that enter our computations, 
in sections 3 and 4 we will drop the family index from ^ without ambiguity. 

The term edY^Y^^ entering the expression for in eq. ([2]) also induces a flavor-changing in- 
teraction with the down quarks for the neutral component of the second Higgs doublet. However, 
this interaction does not affect the computation of T{Z — )• bb), and its contribution to i? — )■ 7 
is negligible with respect to the charged-Higgs contribution as long as edA^/Au <^ (v/mh)'^. Other 
FCNC processes such as BB mixing would put bounds on the combination A^, but this will not be 
relevant to the discussion that follows. 

In the notation of eq. ([3]), the type-I and type-II models are specified by equal to the identity 
matrix in color space, and by the real (and family-universal) coefficients 

Ai = Aii = l/tan/3 (type I), (5) 
4 = -1/4 = l/tan/3 (type II), (6) 

where tan f3 is the ratio of the vevs of the two Higgs doublets in the basis where each of the quark 
mass matrices is obtained from a single Higgs field. 

According to our discussion, the MFV hypothesis includes two other possibilities, namely color- 
singlet and color-octet Higgs doublet that couple to the quarks with arbitrary coefficientil] ^. We 
are going to refer to the first category (singlet) as type-Ill model, while the second possibility (octet) 
will be called type-C model. 

^Even models with generic Yukawa matrices not satisfying the MFV hypothesis show a structure of coupHngs with 
arbitrary coefficients Al^^^. To be phenomenologically viable, the dangerous FCNC effects should be sufficiently suppressed 
via some specific assumption like, e.g., a specific texture of the Yukawa matrices [24j . 
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3 Charged-Higgs contribution to including QCD corrections 



We begin by discussing the radiatively corrected partial decay width of the Z boson in a quark- 
antiquark pair in a model with an additional Higgs doublet. We write it as 



(7) 



where Nc is the color factor {Nc = 3), 5^^. i?) left-handed and right-handed Zqq couplings, 

written in terms of the radiative parameter pq and the radiatively corrected sine of the Weinberg 
angle s^r as (Tg is the third component of the weak isospin, Qq is the electric charge in unit e) 
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Pq - Qq {S 



w> 



(8) 



while the factor Kq contains the QCD, QED and quark- mass corrections. The latter has been com- 
puted up to 0{al) in ref. [25], and at the lowest order it reads: 



3«s I p2 3« 3 Hq 



+ 0(ai 



(9) 



where pq = m^/m^ and Cp = 4/3. 

We assume that the oblique corrections due to the second doublet are negligible, as happens when 
the spectrum of the additional states is approximately custodially symmetric. Then the effect of the 
second Higgs doublet is concentrated in the vertex corrections to T{Z — )• bb). Thus, defining 



SM + H^w)'^ 



(10) 
(11) 



we have 5p(^q^fy-^ = 6{s^f'''^)'^ = 0. In the limit of neglecting the mass of the Z boson with respect to 
the masses of the top quark and the charged Higgs boson we find 
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2 Q, At: si, ^ 



h{th) + ^[CFf2{th) + Clh{th)) 



(13) 



where C]^ = l,C\ = Q [C\ = Cp, C\ = Ca = N^] for Higgs fields in the (1, 2)1/2 [(8, 2)1/2] represen- 
tation, and we omit an overall factor \Vtb\'^ ~ 1. In eqs. (|12p and (|13p th = rh^/m?^^, where rhq is the 
MS quark mass at the scale p and ?7i^+ is the OS mass. The explicit expressions for the functions 
fi{x) are 



X 



re Inx 
T ~ (x- 1)^ 



(14) 
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6x(x-2)^ ( \\ x(-27 + llx) a;(25-9x)lnx 
h\x) = -. T\^Li2 1 H -. — ^ h 



(X-1)2 ^ V Xj (X-1)2 (x-l)3 

6x(3 — x) 12xlnx\ m; 



, 3x ^ / 1\ 3x(l-2x + x2 + ln2x) 6xlnx , , 

where the last term (— 3/i) in the function /2 is introduced to avoid double counting due to the cor- 
rection factor Kq in eq. ([7]), and in the thH~^ coupling we have also kept the contribution proportional 
to the bottom masj§. 

The one-loop terms in 5pb and 5{s\yY agree with the resultjf] of refs. [T71 [22]. The two-loop 
terms were obtained following the lines of the analogous SM calculation that was performed by several 
groups, via different methods, in the early nineties |26j . The SM correction can be actually obtained 
by considering the SM Lagrangian in the limit of vanishing gauge coupling constants, the so-called 
gaugeless limit of the SM [27j. In this limit the gauge bosons play the role of external sources, and 
the propagating fields are those of a Yukawa theory with massless Goldstone bosons. Indeed, in the 
limit mfi+ — 0, the 0{aas'm^ / rn\r) corrections in eq. (fT2]) agree with the known SM result. 

To express the corrections in terms of the OS top mass m^, we must expand the MS top mass 
entering the one-loop part as rht = mt + Smt, with 

6mt = (sin ^ - 4^ mt . (17) 

For the terms proportional to l^nP and \Ad\'^ in eqs. ([T2]) and (fT3]) . this amounts to making the 

and fl. 



substitution rht ~^ n^t and replacing the function /2 with the OS counterparts and /2 , respectively: 



f2{x) = f2{x) 



f f \ , 9/i(x) 
/i(x) + x- 



dx 



N ^ / N Svr 6mt dfi(x) 

/|(x) = /2(x) + —- -x-^. (18) 

OsCp rut ox 



The explicit dependence on the renormalization scale cancels out in the function /g*: 

„„, , 6x(x — 2)^ / 1\ 3x 9xlnx „ , , , , 

whereas has a residual dependence on /x which is compensated for by the implicit scale dependence 
of the ml entering the one- loop parts of eqs. ([12]) and ([13]). Using the OS bottom mass m^, would 
remove this residual scale dependence, but it would introduce large logarithms of the ratio mt/mt in 
the two-loop part of the corrections. 

In fig. [T] we plot, as a function of m^+, the ratio 5^^^ between the two-loop and one-loop con- 
tributions in the terms proportional to \Ay\'^ in eq. (jl2p . for the two cases of color-singlet (lower set 
of lines) and color-octet (upper set of lines) charged Higgs boson. For each case, we show S^'-"^ as 

^Terms proportional to mj,, relevant only for very large values of A^, can also arise from vertices with neutral scalars. 
^In ref. [22] there is a misprint in the overall normalization of the 5g^'^ couplings. 
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Figure 1 : the ratio of two-loop to one- loop charged-Higgs contributions to the Zbb vertex as a function 
of mij+, with the top mass expressed in the MS (solid lines) or OS (dashed lines) renormalization 
scheme. The upper (red) curves are for the model with color-octet Higgs (type C), while the lower 
(blue) curves are for the models with color-singlet Higgs (types I~III). 



obtained using either the central value of the physical (OS) top mass, mt = 173.1 GeV [28j (dashed 
lines), or the corresponding MS value rht{'nit) = 163.5 GeV (solid lines), with the appropriate formulae 
for the two-loop function /2. It can be seen from the figure that, for models of types I-III (i.e. with 
color-singlet charged Higgs) the two- loop corrections are always negative, and they are substantially 
larger when the OS top mass is used in the one-loop part than when the MS mass is used. On the 
other hand, for the model of type C (with color-octet charged Higgs) there is an overall upward shift 
in the two- loop correction due to the additional function /s in eq. (fT2|) . with the result that the sign 
and relative size of the corrections in the OS and MS cases depend on the Higgs mass. For low values 
150 GeV, the two-loop correction approaches zero if the OS top mass is used, and it is positive 
and relatively large if the MS mass is used. Conversely, for larger values ~ 400 GeV the two-loop 
correction approaches zero if the MS top mass is used, and it is negative and relatively large if the 
OS mass is used. As a result, we will see that for 2HDMs of types I-HI a reliable approximation of 
the two-loop result for could be obtained by using the one-loop result expressed in terms of the 
MS top mass. On the other hand, a precise determination of R}, in the 2HDM of type C requires the 
inclusion of the two-loop part of the Zhb vertex correction. 
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From eqs. ([7|)- (|13p we can construct the observable i?;,, which can be written as 
where 

5^ = E ; = U)' + (4)1 (i + Qf£) ■ (21) 

Using the results of ref. [29j we find Sb = 0.6607. Concerning the SM part of Sb, to avoid relying 
indirectly on the measured value of Rf,, we follow ref. [12] and compute it using the values /ff^ = 
0.99426 [29] and (sh/)|m = (1-0063) x 0.23153 = 0.23299, the latter obtained from the measured 
value of sin^ d^^J'j = 0.23153 it 0.00016 corrected for the top-induced contributions specific to the Zbb 
vertejiQ. Finally the factor Cb that includes QCD and mass corrections is obtained from ref. [25j. We 
find, for = 0.118, Cb = 1.0086. 

With the values specified above for the various parameters entering eqs. (|20|) and (|21|) we find a SM 
prediction i^^M ^ 0.21580, nearly la below the measured value Rl""^ = 0.21629 ± 0.00066 ^29]. Since 
the charged-Higgs contributions to g'l and 5^ have the effect of further lowering Rb, stringent bounds 
can be imposed on the parameters mjj+ and by the requirement that the predicted value of Rb in 
a 2HDM be not too far from R'^'^. On the other hand, Rb has little sensitivity on A^, because the 
terms in eqs. (|12p and (jl3p controlled by the latter are suppressed by rhb- Therefore, for the models 
of types III and C in which A^ is a free parameter we will simplify our discussion by setting A^ = 0. 
In the 2HDMs of type I and II the parameter A^ is related to A^ and cannot be set independently 
to zero. However, due to the strong suppression of the contributions controlled by A^, in most of the 
parameter space the predictions of Rb obtained in those two models do not differ significantly from 
the predictions obtained in the type-Ill 2HDM with A^ = 0. More specifically, the predictions of the 
type-I 2HDM, in which A^ = A^, are virtually indistinguishable from those of the type-Ill 2IIDM 
with A^i = for all the values of Au consistent with the measured value of R*^"^. In the 2HDM of type 
II, on the other hand, A^i = —1/Au, and the predictions of Rb differ from the ones obtained in the 
type-Ill 2IIDM with Ad = only for very small values of Ay^. 

Figures [2] and [3] show our determination of Rb in the 2IIDMs of type III and C, respectively, as a 
function of \Au\. In each plot we show two sets of curves for the charged-Higgs mass values = 100 
GeV and ?ti^+ = 400 GeV. In each set, the dashed (solid) curve represents the one-loop (two-loop) 
result expressed in terms of the MS top mass, while the dotted (dot-dashed) curve represents the 
one-loop (two-loop) result expressed in terms of the physical top mass. We also show in each plot the 
measured value R'^'^ (solid horizontal line) and the values la and 2a below (dashed horizontal lines). 
It can be seen that, in both plots, the curves corresponding to the two-loop results (with the top mass 
renormalized either in the MS or in the OS scheme) are practically overlapped. The location of the 
curves corresponding to the one-loop results reflects the behavior that could already be inferred from 



In ref. [12] the correcting factor 1.0063 was not introduced. 
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Color- singlet Higgs 
0.2165 1 < < < < 1 < < < ^ 




Figure 2: Rf^ as a function of \Au\ in the 2HDM with color-singlet Higgs, for = and two different 
values of m^+. The measured value R^^^ = 0.21629 and the values la and 2a below it are displayed 
as horizontal lines. For the meaning of the different curves see the text. 



Color-octet Higgs 
0.2165 1 < < < < 1 < < < 




Figure 3: same as figure [2] in the 2HDM with color-octet Higgs. 
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fig. d} in the 2HDM of type III the one-loop result computed in terms of the MS top mass is very close 
to the two-loop result, while the one- loop result computed in terms of the physical top mass can differ 
significantly. On the other hand, in the 2HDM of type C the quality of the one-loop approximation 
depends on the charged Higgs mass. At low values of , using the physical top mass in the one- loop 
result gives a much better approximation to the two- loop result than using the MS top mass, while the 
situation is reversed at large values of m^+. Therefore, only the use of the two-loop results guarantees 
a precise determination of Rb for all the values of . 

From figures [2] and [3] it is also possible to determine the values of | \ that are disfavored by the 
comparison between R*^"^ and the corresponding theoretical prediction. In the case of the type-Ill 
2HDM, the two-loop curves cross the 2a horizontal line at \Au\ = 0.78 for = 100 GeV and at 
l^nl = 1-35 for mjj+ = 400 GeV. In the case of the type-C 2HDM, they cross it at \Au\ = 0.62 for 
= 100 GeV and at l^^l = 1-10 for = 400 GeV. We checked that the crossing points for 

intermediate values of 171^+ can be determined by linear interpolation of the values given above. 

The predictions of Rb for the type-II 2IIDM (where A^ = —XjA^ = — tan /3) are virtually indistin- 
guishable from those presented in figure [2] for the type-Ill 2IIDM as soon as \A^ > 0.1. The upper 
bounds on \Au\ discussed above translate, for the type-II 2HDM, into \A^ > 1.28 for = 100 GeV 
and \Afi\ > 0.74 for ?ti^+ = 400 GeV. On the other hand, when \Au\ tends to zero the predictions of 
Rb in the type-II 2IIDM decrease quickly, and get two standard deviations below R'^'^ for \Au\ ~ 0.01, 
corresponding to \Ad\ ~ 100. As we will see in the next section, the bounds on Ad coming from the 
process i? — )• Xj 7 can be much stronger than that, but they do not apply to the type-II 2IIDM. 



4 Charged-Higgs contribution to 5 — > 7 at the NLO 

The branching ratio for i? — ?• Xj 7 is fully known at the NLO for a 2IIDM of types I-III. To cover 
also the case of a type-C 2HDM, the only missing ingredient is the determination of the colored-scalar 
0{as) contribution to the Wilson coefficients. We compute it following the analogous computation 
for the type I-II 2IIDM presented in ref. [TU] . 

In the operator basis defined in ref. [IQJ we write the Wilson coefficients at the scale nw, where 
the "full" theory is matched to an effective theory with five quark flavors, as 

= cf\^iw) + 5cf\piw) + W\^iw) + scl^Hm)] , (22) 

47r L J 

where cl''\fj,w) represents the SM contribution {k = 0,1) while 6C^''\iJ,w) represents the charged- 
Higgs contribution. At the LO, the latter is given by 

6Cl'^\fiw) = i = l,...,6, (23) 
6c!f'\m) = Ck{^F^^'{y)-AdA:F:^''{y)] , (24) 

(/i^y) : 






i = 1, 








f\Au\' 


CO 
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where 



(1) _ y{7-^y-Sy^) y\2,y-2) (2) _ yj'i - by) y{3y - 2) 

^y^- 2% -1)3 +4(y-l)4^^^' ^7 ^^^-12(y-l)2 + 6(^^^^^' ^^^^ 

(1) _ y(2 + 5j/-^^) (2) _ ?/(3-?/) y 

^8 ^y^- 8(y-l)3 4(y-l)4^''^' i^i " _ 1)2 2(y-l)3^''^' ^^^^ 

8 - 16(y - 1)2 8(y-l)3^''^' ^8 l^J " _ 1) + _ 1)2 ^^2/' 



with 



y = ^^, (29) 



expressed in terms of the NLO top-quark running mass at the scale fiw and of the OS charged-Higgs 
mass. Again, = 1,C| = [Cj^ = Cf,C\ = Ca = N,] for Higgs fields in the (1, 2)1/2 [(8,2)1/2] 
representation. 

At the NLO, the charged-Higgs contributions to the Wilson coefficients are 

5c\^\m) = z = 1,2,3,5,6, (30) 

6C^^\m) = E'^iy), (31) 



2 

w 

2 

^r^ 
2 



5Ci'\f,w) = Gf{y) + Af{y)ln^ , (32) 



6C^'\^^w) = G^{y) + A^{y)lnJ^ . (33) 



ml. 



The expressions for G^^, A^^g and are rather long and they are reported in the appendix. 
As expected, the dependence in SCj^s cancels out at 0{as) because the functions A^g entering 
eqs. (|32]) and ([33]) satisfy the relation 

where 7q" = 8 is the LO anomalous dimension of the top mass, while ^^j^"^^^ is the matrix of LO 



anomalous dimensions of the Wilson coefficients, whose entries can be found in eq. (8) of ref. 

In fig. m we show the ratio between the NLO and LO charged-Higgs contributions to the Wilson 
coefficients Cy^g as a function of , for both the type-C (color-octet) and type-Ill (color-singlet) 
cases, with the particular choice Ad = (in the color-singlet case this coincides with the type-I 
2HDM). For the type-C 2HDM, the NLO corrections can reach up to ~ 20% of the LO contribution 
at small m^+, and they decrease as the charged-Higgs mass increases, eventually crossing zero. In 
contrast, for the type-Ill 2HDM the two-loop corrections to the Wilson coefficients are always different 
from zero, and have opposite sign with respect to the LO contributions. We checked that the behavior 
of the NLO corrections is qualitatively similar to the one described above even when we allow A^ to 
take on values different from A^. 
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Figure 4: the ratio of NLO to LO charged-Higgs contributions to the Wilson coefficients Cy (solid 
line) and C% (dashed line) as a function of m^+^ with = for the model with color-octet Higgs 
(type C) or color-singlet Higgs (type III). 



The calculation of BR(i? — t- Xs 7) is performed using a modified version of the fortran code SusyBSG 
|30| . The code provides a NLO evaluation of the branching ratio in the MSSM with MFV, including 
the full two-loop gluino contributions to the Wilson coefficients [31j. The current public version (1.3) 
includes also the options of evaluating the branching ratio in the SM, in the type-II 2HDM and in 
the MSSM with two-loop gluino contributions computed in the effective Lagrangian approximation. 
We enlarged the 2HDM option to include also the type-I, type-Ill and type-C models, thus covering 
all four types of 2HDM compatible with MFvi The relation between the Wilson coefficients and 
BR(i3 —7- Xg^) is computed at NLO along the lines of ref. [32], but the free renormalization scales 
entering the NLO calculation are adjusted in such a way as to mimic the Next-to-Next-to-Leading 
Order (NNLO) contributions presented in ref. ^3]. When the SM input parameters are set to the 
partially outdated values used in ref. [33], SusyBSG gives a SM prediction for BR(i? — )• 7) of 
3.15 X 10"'', in full agreement with the NNLO result of that paper. Very good numerical agreement 
is also found with the results of the partial NNLO implementation of the type-II 2HDM in ref. [33], 
which combines NNLO anomalous dimensions and matrix elements with NLO Wilson coefficients. We 
public version of SusyBSG with this new feature will be released soon. 
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A =0.3 A =-0.3 

U LI 




Figure 5: BR(i? — t- 7) as a function of A^i in the type-III (solid line) and type-C (dashed line) 
models for Au = 0.3 (left panel) and Au = —0.3 (right panel) and two different values of mjj+. The 
horizontal dashed lines specify the 2a interval around the experimental value of BR(i3 — t- Xsj). 



take into account a recent update [34] in the calculation of the normalization factor for the branching 
ratio as well as the latest central value of the top mass [28], which results in a modest enhancement 
of the SM prediction for BR{B X., 7) to 3.28 x 10"^. 

In the 2HDMs of types I and II, the requirement of consistency between the theoretical prediction 
and the measured value of BR{B — )• Xs 7) allows us to set bounds on the parameters and tan /3, 
the latter determining both Higgs-quark couplings Au and A^. More specifically, in the type-I 2IIDM 
the charged-Higgs contribution to the Wilson coefficients Cy^g scales like 1 / tan , therefore it is 
possible to derive, for each given value of m^+, a lower bound on tan/3 (i.e., an upper bound on 

= Ad) which is however much less stringent than the corresponding bound derived from Rf^. In 
the type-II 2HDM there is a tan /3-independent contribution to the Wilson coefficient, which allows to 
set an absolute lower bound on m^+. These bounds have been extensively discussed in the literature 
(see, e.g., refs. [11^ [121 113j ) and we will not further consider them here. In the models of type III 
and C, on the other hand, the parameters Au and A^ are unrelated to each other. As can be seen in 
eqs. dMl) and (p5|) . the charged-Higgs contributions to the Wilson coefficients 6*7^8 include two terms 
controlled by \Au\^ and A^A^^, respectively. The bounds on A^ derived in the previous section tell us 
that the former term cannot be too large, while the latter can be significant for large values of A^. 
Furthermore, its effect on the branching ratio depends on the relative sign between A^ and A^. 

In fig. [5] we show BR(i? — )• Xg 7) as a function of A^-, for both the type-III (solid lines) and type-C 
(dashed lines) 2IIDM, and for the representative choices A^ = ±0.3 (the latter are allowed by R^, as 
can be seen in figs. [2] and [3|). The left panel displays the case of same sign between Au and A^, while 
the right panel shows the case of opposite sign. Each plot contains two sets of curves corresponding 
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A =0.3 



A 

pq 



— Type III 
-- TypeC 




Figure 6: BR(5 
models for Au = 
interval around the experimental value of BR(i? — )• Xj 7). 



— )• 7) as a function of in the type-Ill (solid line) and type-C (dashed line) 
0.3 and two different values of m^+. The horizontal dashed lines specify the 2a 



to = 100 GeV and = 400 GeV, respectively. The horizontal dashed lines mark the 95% 

C.L. band around the experimental value BR(5 ^ X, 7) = (3.52 ± 0.25) x 10""^ [35j. The band 
also includes, added in quadrature, the theoretical error on the 2HDM prediction (we conservatively 
estimate this error as 10% of the SM prediction for the branching ratio). From fig. [5] it is clear that, 
unless \Au\ is extremely small, the process B ^ Xg'j sets stringent limits on A^i. Focusing on the 
case of color-singlet Higgs we see that, for A^ = 0.3, the values of A^ that allow for a branching ratio 
inside the 95% C.L. band are Ad < 0.9 for = 100 GeV and Ad < 2.5 for = 400 GeV; for 
Au = —0.3 the bounds are a little less stringent, i.e. Ad < 1.3 for mjj+ = 100 GeV and Ad < 4 for 
= 400 GeV. From the figure it is also apparent that the bound on Ad for a fixed value of Au 
is almost independent of the colored or colorless nature of the charged Higgs, with the colored case 
showing only slightly stronger bounds. However, as seen in the previous section, the bounds on Au 
derived from Rb are more dependent on the nature of the Higgs, so that the allowed regions for the 
Ad coefficient are in fact different for the color-singlet and color-octet charged Higgs. 

The case of same sign between A^ and Ad has the peculiarity that, as shown in fig. [6l there are 
actually two ranges of values of Ad that fit inside the 2a allowed band for BR(i? — )• Xg'j). This 
is related to the fact that, in this case, the sign of the charged-Higgs contribution to the Wilson 
coefficient of the magnetic dipole operator, , is opposite to the sign of the SM contribution, Cj^. 
Since BR(i? — t- Xg 7) is roughly proportional to \Cj^ + C^^p, as Ad increases the branching ratio goes 
to zero when ~ — Cy^, and it goes back inside the 2(T-allowed band when ~ —20^^. Thus, 
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the two ranges of possible values of differ by the sign of the amplitude A(b — t- 57), basically the sign 
of the Wilson coefficient C7. Although i? — ?• X, 7 allows both ranges of values for A^i, there are other 
observables that are sensitive to the sign of C7, thus selecting one of the two options. Among them, 
we cite BR(i? — t- Xsl~^l~) [36] and the isospin-breaking asymmetry that can be constructed from 
the exclusive neutral and charged B — K*^ decay modes [37l [16]. Although, for both observables, 
neither the experimental result nor the theoretical prediction is at the same level of accuracy as for 
BR(S — )• -^^7), these observables still give a compelling indication that the sign of A{b — )• 57) is that 
of the SM contribution, thus eliminating the large- option. 

5 Conclusions 

Minimal flavor violation is a very popular criterion that is used to suppress FCNC effects in models 
with new particles at the TeV scale. The enforcement of the MFV criterion to the simplest extension 
of the SM, i.e. a model with a second Higgs doublet, allows the possibility of color-singlet or color-octet 
Higgs field. For both cases we have considered the two-loop QCD corrections to the charged Higgs 
boson contribution to Rh. We found that for all four types of 2HDM with MFV, after the inclusion 
of the two-loop QCD corrections, the prediction for this observable is practically insensitive to the 
the choice of renormalization scheme for the top mass entering the one-loop part of the calculation. 
Thus, the upper bound on the coupling \Au\ derived from Rf, is improved, which for the type-I and 
type-II models translates in a more robust lower bound on tan/3. We have also computed the 0{as) 
contributions to the Wilson coefficients relevant to the process B ^ Xg'j for the color-octet case. 
This was the last missing ingredient to obtain a determination of BR(i? — )• 7) at the NLO level 
for all 2HDM with MFV. After the inclusion of the NLO corrections it is found that, in the region 
allowed by the present experimental results, the i? — t- X, 7 transition is fairly insensitive to the colored 
or colorless nature of the charged Higgs. Furthermore, in type-Ill and C models, the bounds on the 
Ad parameter that can be obtained from BR(i? — >• Xg 7) and other observables rule out the large- 
region, where effects proportional to the bottom mass could become important. 
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Appendix: Analytical expressions for the NLO Wilson coefficients 

In this appendix we report the analytic expressions for the functions Gy^g and A^^g entering SCj^g. 
We find 



^"-'^ V - y) + 3(,-i)4 1- y 



+ \Au\'^y 



2(-3-j/ + 12y^-2y3) 
3(y-l)4 

y(18 - 37y + 8y^ 



Iny + 



1 



+ 
■i 



6(y-l)4 yj ' 6(y-l)5 

50 + 251y - 174?/2 _ I92y3 + 21y^ 
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, 2/(-14 + 23y + 3y2) ^ 



In y 
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797 - 5436y + 756%^ - 1202y3 



Iny 



+ CkCl^AdAly 



648(y - 1)4 
-l^ + ^^^.LiJl-i,+ 



[I8(y-1)^ 
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18(y-l)3 
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(Al) 
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12(y-l)2 6(y-l) 
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+ \Aur-^y 
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(A5) 



The results for type I-III models are recovered setting = 1 and = 0, while the result for 
the case of a colored charged scalar in the adjoint of SU{3) is obtained with = Cp and = N^,. 



Note added: after the publication of our paper, we became aware of ref. [38], which contains 
two-loop formulae for the Wilson coefficients of the operators relevant to 6 — )• 57 in a generic model with 
a heavy fermion and a heavy scalar. While ref. [38] does not specifically discuss the case of the type-C 
2HDM, the results presented in our appendix can be reproduced with appropriate substitutions in the 
results of ref. [38] , taking into account the different renormalization conditions on the charged-Higgs 
mass. We thank Christoph Bobeth for performing this useful comparison. 
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